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Abstract. We present the results of a high angular res- 
olution, multi-transition analysis of the molecular gas in 
M82. The analysis is based on the two lowest transitions 
of 12 CO and the ground transition of the rare isotopes 
13 CO and C 18 measured with the PdBI, the BIMA ar- 
ray and the IRAM 30m telescope. 

In order to address the question of how the intrinsic molec- 
ular cloud properties are influenced by massive star for- 
mation we have carried out radiative transfer calculations 
based on the observed CO line ratios. The calculations 
suggest that the kinetic temperature of the molecular gas 
is high in regions with strong star formation and drops 
towards the outer molecular lobes with less ongoing star 
formation. The location of the highest kinetic temperature 
is coincident with that of the mid infrared (MIR) peaks 
which trace emission from hot dust. The hot gas is associ- 
ated with low H2 densities while the cold gas in the outer 
molecular lobes has high H2 densities. We find that CO 
intensities do not trace H2 column densities well. Most of 
the molecular gas is distributed in a double-lobed distri- 
bution which surrounds the starburst. A detailed analysis 
of the conversion factor from CO intensity to H2 column 
density shows that Xco depends on the excitation con- 
ditions. We find Xco ~ T^ 1 n(H 2 ) 1 ^ 2 , as expected for 
virialized clouds. 
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1. Introduction 

M82 is rega rded as the archetypical starburst galaxy 
(Rieke etal. 1980 ). Its distance of only 3.9 Mpc (Sakai 
& Madore 1999|) makes M82 an excellent laboratory for 
studying the relevant physical processes connected with 
starburst activity in detail. The central few hundred par- 
sees of this galaxy are heavily obscured by dust and 
gas which hides the central starburst region against di- 



rect observations at optical wavelengths. Evidence for the 
strong star-forming activity in the central region comes 



from radio (e.g. Kronberg etal. 1985; Wills etal 



1999) 



and infrared observations (e.g. Telesco & Gezari |1992|) 
and also from the prominent bipolar outflow visible in 
Ha (e.g. Bland & Tully |l988[ McKeith et al. |l995| ; Shop- 
bell fc Bla nd-Hawthorn 1998 ) and X-rays (e.g. Bregman 
etal. |!995p . The massive star formation (SF) is believed 
to be fueled by the large amount of molecular gas which 
is present in the centre of M82. 

On the other hand, SF affects the distribution, kinemat- 
ics and physical conditions of the surrounding interstel- 
lar medium (ISM). Early studies of the distribution of 
the molecular gas in M82 unveiled a double-lobed cir- 
cumnuclear distribution of CO which was interpreted as 
a molecular torus with a depletion of molecular gas in 
the central region (Nakai etal. 1987). More recent high- 
resolution studies by Shen & Lo (1995) using BIMA and 
Neininger etal. ( 1998| ) using the IRAM interferometer at 
Plateau de Bure (PdBI) showed a third molecular peak 
65 pc west of the 2.2/im nucleus (Dietz etal. 1986). Using 



these high-resolution CO maps WeiB etal. (1999) iden- 



tified an expanding superbubble in the molecular gas of 
M 82 which linked the triple peak CO distribution and its 
disturbed velocity field to the prominent outflow visible in 
Ha and X-rays. 

Multi-transition analyses of molecular emission lines (CO, 
CS, HCN) showed that the starburst also affects the phys- 
ical conditions of th e molecular gas (W ild et al. |1992| 



Hcnkel & Bally |1985| ; Brouillet & Schilke [1993J) . A large 
fraction of the molecular gas is concentrated in warm 
(T km = 50 K) and dense ( n(H 2 ) = 10 4 cm" 3 ) clouds (Wild 
etal. 1995; Gusten etal. 1993). In a recent study Mao 
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et al. (200C) analyzed CO mm and sub-mm emission lines 
in M 82 to investigate the physical properties of the molec- 
ular clouds. They conclude that the bulk of CO emission 
arises from photon-dominated regions (PDRs) while trac- 
ers of high-density gas like CS and HCN are less affected 
by the strong UV radiation from massive stars. 
Even though these studies already provided a good global 
picture of the physical conditions of the molecular gas in 
M82, no detailed high-resolution study exists so far that 
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allows to resolve variations of the excitation conditions of 
the molecular gas over the central part of M82. In this 
paper we present the results of a high angular resolution, 
multi-transition CO analysis and compare the intrinsic 
gas properties with observations of high - level star forma- 
tion. In Sect. 2 we briefly summarize our observations and 
the data reduction. In Sect. 3 we describe the main re- 
sults including a description of the CO morphology and 
kinematics, the CO line ratios, results from the LVG cal- 
culation and on Xqo- I n Sect. 4 we compare our results to 
previous studies. Sect. 5 summarizes our conclusions. 

2. Observations 

2.1. PdBI Observations 

We have used the PdBI to observe the 12 CO(J = 2 1) 
(v obs = 230.361GHz) and C 18 0(J = 1^0) 
{v bs — 109.698 GHz) emission lines in the central 
region of M 82. The observations were carried out in April 
1997. Due to the dual frequency setup of the PdBI we 
were able to observe both emission lines simultaneously. 
The observations were carried out in mosaic mode with 
seven pointings covering the central kpc of M82. The 
central pointing was centred on the 2.2/im nucleus at 
a = 09^55 m 51 s .94,<5 = 69°40'47.14" (J2000.0) (Dietz 
etal. 1986). The other pointings were shifted with 



respect to the central position by (a, 5) = (— 30", — 8"), 
(-20", -6"), (-10", -4"), (10", 4"), (20", 8"), (30", 10") 
which ensured sufficient overlap of the observed fields at 
230 GHz. The primary beam of the PdBI is 22" and 45" 
at 230 GHz and 109 GHz respectively. The observations 
were carried out in the DC2 antenna configuration with 
baselines ranging from 24 m to 176 m leading to a synthe- 
sized beam of 3.8" x 3.5" at 109 GHz and 1.5" x 1.4" at 
230 GHz. 

The C 18 0(J= 1 — > 0) data were recorded using two cor- 
relator units leading to a total bandwidth of 780 kms -1 
with 6.83 kms" 1 resolution. For the 12 CO(J = l-»0) 
transition we used four correlator units which resulted in 
a total bandwidth of 390 kms -1 and a velocity resolution 
of 3.25 kms" 1 . The 12 CO(J = 2 -► 1) emission line 
was observed in the lower sideband of the 230 GHz, the 
C 18 0(J= 1 — > 0) emission line in the upper sideband of 
the 109 GHz receiver. The flux and complex bandpass 
calibration was determined by observing 3C273 and 
MWC349. The nearby calibrator 0836+710 was used 
as a secondary amplitude and phase calibrator. The 
seven fields were combined in a mosaic and subsequently 
cleaned using the MAPPING procedure of the GILDAS 
software package. This yields a roughly constant sensi- 
tivity along the major axis of M 82 with an rms noise of 
6mJy/beam at 109 GHz and 30mJy/beam at 230 GHz. 
For both data sets the channels with vi sr > 385 kms -1 
and vi sr < 30 kms -1 were used to generate a continuum 
map at 109 GHz and 230 GHz. The continuum emission 
was subtracted from both emission line data cubes. 



2.2. IRAM 30m Telescope Observations 

In addition to the high-resolution CO data we observed 
the 12 CO( J = 1 -» 0), 12 CO(J = 2^1) and 13 CO( J = 

1 — > 0) emission lines with the IRAM 30m telescope in on- 
the-fly mode. The observations covered an area of 3' x 3' 
centred on the 2.2^m nucleus. The 12 CO(J=l — > 0) ob- 
servations were carried out in Nov 1997. The 12 CO(J = 

2 — > 1) and 13 CO(J = 1 — * 0) data were observed in 
Nov. 1997, Dec. 1998 and June 1999. For all observations 
we used the same observing procedure: The scanning ve- 
locity was 2"/sec and the readout sampling lsec lead- 
ing to a spatial separation of 2" between individual spec- 
tra in scanning direction. The spatial separation between 
individual scans was 4". Thus each on-the-fly map was 
sampled on a 2" x 4" grid. For the 12 CO(J = 1 -> 0) 
transition we performed two coverages, for the other two 
transitions we performed four coverages with perpendicu- 
lar scanning directions. The combined data therefore were 
sampled on a 2" x 2" grid. After first-order baseline sub- 
traction the spectra were summed on a 3" x 3" grid us- 
ing the beam (13" at 230 GHz, 22" at 115 GHz) and the 
rms noise level for weighting. This observing and reduc- 
tion procedure guarantees a smooth data sampling and 
avoids artifacts in the combination with the interferomet- 
ric data. The total integration time per beam was 65 sec 
for the 12 CO(J=l -> 0), 45 sec for the 12 CO(J = 2 -► 1) 
and 130 sec for 13 CO( J = 1 — * 0) transition resulting in 
an rms noise of 40 mK, 65 mK and 15 mK. As backends 
we used the autocorrelators which lead to a total band- 
width and velocity resolution of 650 kms /2.6 kms" 1 , 
650 kms _1 /1.3 kms _1 and 695 kms~ 1 /2.7 kms _1 for the 
12 CO(J = l -* 0), 12 CO(J = 2^ 1) and 13 CO(J = l -» 0) 
transitions respectively. For the conversion from Tjj[ t o T m b 
we used F eff /B eff = 1.35 at 115 GHz (Guelin e tal. |l995|) 
and Feff/Beff = 2.05 at 230 GHz (Greve etal. |l998|) . 



2.3. Short-Spacing Correction 

To ensure that the interferometric line intensities do not 
suffer from missing flux due to extended emission we com- 
bined the interferometer and the single-dish data cubes. 
For the combination we used a method that works on the 
final reduced (CLEANed and corrected for primary beam 
attenuation) interferometer cubes. The only free param- 
eter in this method is the choice of which part of the 
Mv-plane in the interferometer cube is replaced by the 
single-dish values. A detailed description of the method is 
given in appendix A. The parameters for the 30m beam 
sizes, the corresponding effective diameter of the 30m tele- 
scope, the shortest baseline, the replaced part of the uv- 
plane and the missing flux of the interferometer maps are 
given in Table [l]. All reduction steps were done using the 
MIRIAD software package. We applied the short-spacing 
correction to the 12 CO( J = 1 — > 0) cube obtained by Shen 
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Table 1. Summary of the relevant parameters of the 
short-spacing correction. 



2 CO (1-0) 



2 CO(2-l) 



'CO (1-0) 



FWHM 

Deff 
Smin 

SD 



D, 



k\ 

miss, flux 



22" 
28.1m 
unknown 
<9.4 
20 % 



11" 
24.4 m 
24 m 
<18.5 
60 % 



22" 
28.1m 
24 m 
<8.0 

35 % 



From top to bottom the parameters are the FWHM of the 
IRAM 30 m telescope beam, the corresponding effective area 
of the 30 m telescope, the shortest baseline in the interferome- 
ter data, the part of the airplane that has been replaced by the 
single-dish data and the average missing flux in the interferom- 
eter maps. Since the shortest baseline of the 12 CO(J=l — > 0) 
observations from Shen & Lo is not given in their paper we 
replaced the part of the tti>-plane which corresponds to the 
effective diameter of the 30 m telescope. 



& Lo pgga ), the 13 CO(J=l -> 
etal. (|99§ and to the 12 CO(J = 



0) cube from Neininger 
= 2 -> 1) cube. 



3. Results 



3.1. CO morphology and kinematics 

Figs. and | show the integrated 12 CO( J = 2 -> 1) and 
C 18 0(J=1 — + 0) line intensities. The overall morphology 
of both images is very similar to the 12 CO( J = 1 — > 0) 
distribution published by Shen & Lo (1995) and the 
13 CO(J = 1 — > 0) distribution published by Neininger 
etal. (199S). It shows a triple peak morphology of which 



the two outer lobes have been interpreted as the edge of 
a cen tral molecular toroid (Nakai et al. 1987] , Shen & Lo 



1995 ) and a weaker central peak located 65 pc west of the 
M82's centre (2.2/im peak; Dietz etal. |l986j ). The two 
outer lobes have a projected separation of 410 pc (26"). 
The separation of the central and the western molecular 
lobe is only about 130 pc (8"). More diffuse CO emission 
is detected in the 12 CO(J = 2 — > 1) intensity distribu- 
tion east and west of the CO peaks and in the south-west 
of the galaxy. The eastern part of the CO distribution is 
significantly warped to the north. The total extent of the 
emission region is about 1 kpc from east to west. With re- 
spect to M 82 's centre the distribution of the molecular gas 
is clearly displaced to the west. South of the central and 
western CO peak two CO spurs are detected (see Fig|l|). 
They extend about 100 pc below the main molecular disk 
and join just below the expanding molecular superbubble 
which is located between the central and western CO peak 
(Neininger etal. [1998| Weifi etal. |1999| ). At the same lo- 
cation hot gas emerges into the halo of M 82 (e.g. Shopbell 
& Bland-Hawthorn 199S; Bregman etal. 1995) support- 



ing the idea that the CO spurs indicate the walls of the 
superbubble. 

Note that the chain of CO emission south of the eastern 
end of the 12 CO(J = 2 — > 1) distribution is most likely not 




Fig. 3. A pv-diagram along the major axis of M 82 in the 
12 CO(J = 2 — » 1) transition. The contours correspond to 
4, 8, 12, 14, 18, and 22 K. The pv-diagram is centred on 
the brightest radio continuum point source SNR 41.9+58 
(see Figs, [j] & ||). The velocity crowded regimes at 20", 
5" and —7" offset correspond to the western, central and 
eastern CO peak. 



real but an artifact from the primary beam correction. 
The kinematic of the central 400 pc is dominated by 
solid body rotation. The rotation amplitude is about 
200 kms -1 ranging from 115 kms -1 at the western peak 
up to 320 km s _1 at the eastern peak. A pv-diagram along 
the major axis of M82 in the 12 CO(J = 2 — » 1) transition 
is shown in Fig. |[ (For the corresponding diagram in the 
C 18 0(J = 1 0) data see Weifi etal. [L999D . The pv- 
diagram is centred on the brightest supernova remnant 
SNR 41.9+58. The intense, velocity crowded regions at 
20", 5" and —7" offset correspond to the western, cen- 
tral and eastern CO peak. Between the central and west- 
ern CO peaks two velocity components at 100 kms -1 and 
190 kms -1 are detected. These features have been inter- 
preted as an expanding superbubble. The velocity of the 
CO spurs is about 140 kms -1 (see Figs. ^ and ||) which 
is similar to the centroid velocity of the expanding su- 
perbubble. Outside the central 400 pc the CO rotation 
curve flattens. The dynamical centre derived from the 
12 CO(J = 2^1) and C 18 0(J = 1^0) data agrees 
very well with the value of v sys = 225 ± 10 kms -1 pub- 
lished by Shen & L o (|1995|) for the 12 CO (J = 1 -> 0), 
and Neininger etal. fll998| ) for the 13 CO(J=l -> 0) tran- 
sition. The channclmaps of the 12 CO(J = 2 — ► 1) and 
C 18 0( J = 1 — > 0) line emission are presented in Figs. |^ 
and ||. 



3.2. Location of the starburst with respect to CO 

Most tracers of star formation in M 82 indicate that the 
highest star-forming activity is not associated with the 
molecular peaks, which presumably indicate the location 
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9 h 55 m 55 s 9 h 55 m 50 s 
RA J2000 



Fig. 1. Integrated 12 CO(J = 2 — > 1) line intensity derived from the short-spacing corrected data cube. The contours 
correspond to 200, 400, 600, 800, 1200, 1600, 2000 and 2400 Kkms" 1 . The cross indicates the centre of M82 (2.2^m 
peak), the star represents the location of the brightest radio continuum point source (SNR 41.9+58) in M82. The 
dotted line shows the 50% sensitivity level of the primary beam mosaic. 



69°41'00" - 



o 69°40'50" - 
—) 

o 

LJ 

Q 

69°40'40" - 



69 o 40'30" - 




RA J2000 



Fig. 2. Integrated C 18 0(J=1 -> 0) line intensity. The contours correspond to 2, 8, 14, 20, 26 and 32 Kkrns" 1 . The 
cross indicates the centre of M82 (2.2/im peak), the star marks the location of the brightest radio continuum point 
source (SNR 41.9+58) in M 82. 



of the reservoirs for the 'fuel' for star formation, but rather 
takes place between the peaks. The high-resolution 12.4 
/um image of the central region of M 82 published by Tele- 
sco & Gezari ( |1992| ) suggests that the young stellar clus- 
ters, which heat the dust, are located between the west- 
ern molecular lobe and the 2.2 /xm nucleus (western mid 
infrared (MIR) peaks), at the central CO peak, and be- 
tween the central CO peak and the eastern CO lobe (east- 
ern MIR peak). A similar morphology is visible in the 



Nell line emission (Achtermann & Lacy 1995). The radio 



continuum point sources, which are believed to be super- 
nova remnants (SNR) and compact HII regions, are spread 



across a much wider region and seem to avoid MIR and 
Nell peaks (Kronberg et al. |1985| ). Only the strongest SNR 
in M 82, SNR41.9+58, appears to be related to features at 
other wavelengths: it is located near the centre of the ex- 
panding molecular superbubble, between the central and 
western CO peak, from which hot X-ray emitting gas is re- 



leased into the halo of M 82 (Weifi et al. |1999| ). At the same 
location recent radio continuum studies by Wills et al. 
( 1999| ) identified a blow-out in the form of a cone of miss- 
ing 5-GHz continuum emission. In the same study three 
other chimneys were identified within the central 300 pc 
of M 82. All these observations indicate that the regions of 



A. Weifi etal.: The state of the molecular gas in M82 5 



20 
10 


-10 
-20 



1 1 I 1 1 1 I 



T 1 1 1 1 1 1 - 




I I I I I I I I I I I I I I I 

250.6 



I I I — | I I I | I I I | I I I 



SO 1 



I I I | I I I | I I I | I I I 




++++ 



-EE*- 

I I I | I I I | I I I | I I 



+ 



-i — i — i — i— 



-i — i — i — i — r- 



+ 



1 1 1 — 1 1 1 1 1 1 — 1 1 1 1 1 1 




1 1 1 — 1 1 1 1 1 1 — 1 1 1 1 1 1 

„ 237.8 



I I I — | I I I | I — I I | I I I 



eg" 



I I I | i i i | i i i | i i i 




++++ 



i i | i i i | i i i | i i i 



0,0 



:: '. - -'^S^>'* 



+ 



I I I | I I I | I — I I | I I I 




' ' ' I ' ' ' I ' ' ' I ' ' ' 
225 



i i i | i i i | i — i i | i i i 



i i i | i i i | i i i | i i i 



Mi l l 




I I I | I I I | I I I | I I I 



i i i I i i i I i — i i | i i i 




i i i I i i i I i — i i | i i i 



i i i I i i i I i — i i | i i i 




i i i I i i i I i i i I i i i 



Mill 




I 



_l I I 1_ 



+ 



20.20 [K] 
20 



10 







40 20 -20 -40 

Fig. 4. 12 CO( J = 2 — * 1) velocity channel maps of the short-spacing corrected data cube. Offsets in RA, DEC are given 
relative to the 2.2y«m nucleus indicated by the cross in each channel map. The velocity spacing between individual 
maps is 12.8 kms -1 , the central velocity (Vi sr ) of each map is given in the top right corner in units of kms . The 
intensities are corrected for attenuation of the primary beam. The contours correspond to 6, 12, 18, 24 and 30 K. The 
dotted line in each channel map shows the 50% sensitivity level from the primary beam mosaic. 



violent star formation are confined by the molecular lobes. 
Since no indications for high activity have been found at 
the 2.2 fim nucleus itself, it seems that the star burst is 
arranged in a toroidal topology around the nucleus. 

3.3. CO line ratios 

To calculate the line ratios properly we used the short- 
spacing corrected 12 CO(J=l -► 0), 12 CO(J = 2 -> 1) and 
13 CO( J — 1 — > 0) data cubes. Note that the missing flux in 
the pure interferometric maps can be as high as 60% (see 
Table . Therefore the short-spacing correction is vital to 
derive proper line ratios. The short-spacing correction is 
less crucial for the peak intensities. Here the missing flux 



is 10%-30% only. The 12 CO(J=l -» 0), 12 CO(J = 2 -» 1) 
and C 18 0(J= 1 — > 0) data were smoothed to the resolu- 
tion of the 13 CO(J=l -> 0) observations (4.2"). Since no 
single-dish data were obtained for the C ls O(J = 1 — > 0) 
transition we applied the missing flux factors derived from 
the 13 CO( J — 1 — > 0) peak intensity distribution to the 
C ls O(J = 1 — > 0) observations. This procedure is justi- 
fied, because the frequency of both transitions is similar 
and the observations were carried out in the same configu- 
rations with the PdBI. This leads to similar vy-coverages 
for both observations. Furthermore the morphology in the 
interferometer maps is similar and both transitions are 
optically thin (see Sect. 3.4). To take the remaining un- 
certainties into account we assumed an error of 50% for 
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Fig. 5. C 0(J=1 — * 0) velocity channel maps of the PdBI data cube. The figure has the same layout and specifica- 
tions as Fig. [| Contours correspond to 0.05, 0.15, 0.25, 0.35, 0.45 and 0.55 K. 



the C 18 0(J=1 — + 0) line intensities. The line ratios were 
calculated at 19 positions across the CO distribution of 
M82. The spacing between individual positions is about 
4". The analyzed positions are marked by the crosses in 
Fig. | The circles indicate the FWHM of 4.2" used in 
the study. The positions include all molecular peaks, the 
2.2/im nucleus, the MIR peaks, the CO spurs and the 
diffuse emission in the outer regions of M82. For clarity 
the positions have been labeled 1 to 16 from east to west. 
Positions 17 to 19 correspond to positions on the CO out- 
flow (see Fig. ||). The line ratios at the analyzed positions 
are summarized in Tab Errors include 10% uncertainty 
of the flux calibrators, errors of the amplitude calibration 
(typically about 10%) and statistical errors. Our high- 
resolution line ratios for 12 CO and 13 CO differ slightly 
from values derived from single dish observations by Mao 
etal. (|2000D. But our data confirms that 12 CO(J = 2 -> 



l)/ 12 CO(J = 1 -> 0) ratios lar ger th an 1.8 (e.g. Knapp 
et al. |1980[ Olofsson & Rydbeck [1984 Loiseau et al. [1990| ) 
can firmly be rejected. 12 CO(J=l -> 0)/ 13 CO(J=1 -► 0) 
and 12 CO(J=l -> 0)/C 18 O(J=l -> 0) line intensity ra- 
tios are about 10-20 and 40-60 respectively. 

3.4- Radiative transfer calculations 

The excitation conditions of the CO-emitting volume were 
modeled using a spherical, isothermal one-component 
large velocity gradient (LVG) model (Goldreich & Kwan 
1974, dc Jong etal. 1975| ). LVG line intensities were cal- 
culated for a kinetic temperature and H2 density range 
from 5 K to 200 K by 5 K and log n(ff 2 ) from 1.8 to 5.0 
by 0.2 respectively. In addition, we varied the CO abun- 
dance relative to H2, [CO], per velocity gradient and the 
fractional 13 CO and C 18 abundances ([CO]/grad(V): 
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Table 2. CO line ratios in M82 at 4" resolution. The offsets are given relative to a = 09' l 55 m 51 s .94, S = 69°40'47.14" 
(J2000.0). For Column 8 & 9 we have corrected the C 18 intensities with the missing flux factor determined from the 
13 CO data. The corresponding percentage is given in Column 10. Errors include 10% uncertainty of the flux calibrators, 
errors of the amplitude calibration (typically about 10%) and statistical errors. For the C 18 intensities we assumed 
an error of 50% due to the unknown missing flux. 



I(CO(2-l)) T(CO(2-l)) I(CO(1-0)) T(CO(1-0)) I(CO(1-0)) T(CO(1-0)) MP^m 1 

I(CO(1-0)) T(CO(1-0)) I( 1: iCO(l-0)) T( la CO(l-0)) I(C 18 O(l-0)) T(C ls O(l-0)) 

[%] 



1 


16.5 


7.5 


1.16 ± 0.3 


1.02 ± 0.3 


15.2 ± 4.9 


19.2 ± 5.7 


44 ± 24 


66 ± 35 


23.5 


2 


14.5 


5.0 


1.17 ± 0.3 


0.96 ± 0.3 


12.7 ± 3.9 


15.4 ± 4.9 


38 ± 21 


45 ± 24 


10.3 


3 


11.5 


3.0 


1.12 ± 0.3 


0.98 ± 0.3 


11.5 ± 3.4 


12.2 ± 3.7 


40 ± 21 


44 ± 24 


8.0 


4 


9.5 


1.0 


1.19 ±0.3 


1.09 ±0.3 


12.0 ±3.6 


13.0 ±4.1 


39 ±21 


43 ±23 


7.0 


5 


6.5 


0.5 


1.34 ±0.4 


1.03 ±0.3 


14.2 ±4.4 


12.8 ±4.2 


37 ±20 


40 ±22 


9.0 


6 


4.0 


0.5 


1.42 ±0.4 


1.28 ±0.4 


15.0 ±4.8 


12.6 ±4.1 


33 ± 18 


37 ±20 


12.4 


7 


2.0 


-0.5 


1.39 ± 0.4 


1.26 ±0.4 


17.3 ±5.7 


16.5 ± 6.1 


37 ±20 


61 ±33 


10.3 


8 


-1.0 


-1.5 


1.06 ±0.3 


1.01 ±0.3 


25.9 ±8.5 


20.7 ±8.0 


36 ±20 


56 ±30 


13.0 


9 


-4.0 


-2.0 


0.98 ±0.3 


0.98 ±0.3 


21.3 ± 6.6 


22.1 ± 8.4 


53 ±28 


56 ±30 


15.6 


10 


-6.5 


-3.0 


1.07 ±0.3 


0.95 ±0.3 


25.4 ±8.4 


17.0 ±6.0 


41 ± 22 


45 ± 24 


12.9 


11 


-10.0 


-4.0 


1.12 ±0.3 


1.08 ±0.3 


13.2 ± 3.9 


11.7 ±3.8 


45 ± 24 


29 ± 15 


14.0 


12 


-14.0 


-4.0 


1.00 ±0.3 


0.98 ±0.3 


15.3 ±4.5 


11.5 ±3.5 


48 ±26 


36 ±20 


4.3 


13 


-17.5 


-4.0 


1.14 ±0.3 


1.02 ±0.3 


14.6 ±4.4 


10.9 ±3.3 


48 ±26 


39 ±21 


6.6 


14 


-20.5 


-4.5 


1.05 ±0.3 


0.95 ±0.3 


14.3 ±4.3 


11.1 ± 3.6 


49 ±27 


39 ±21 


4.1 


15 


-23.5 


-4.5 


1.07 ±0.3 


1.00 ±0.3 


16.3 ± 5.1 


14.7 ±4.7 


60 ±32 


48 ±26 


5.6 


16 


-26.5 


-4.0 


1.14 ±0.3 


1.09 ±0.3 


14.5 ±4.8 


14.3 ±4.9 


62 ± 34 


56 ±30 


6.2 


17 


-2.5 


-5.5 


1.02 ±0.3 


1.07 ±0.3 


20.2 ± 6.8 


21.2 ± 8.3 


40 ±22 


54 ±29 


19.0 


18 


-3.5 


-9.0 


1.20 ±0.3 


0.91 ±0.3 


12.0 ±4.1 


21.5 ± 8.6 


29 ± 16 


52 ± 28 


27.4 


19 


-8.5 


-7.0 


1.36 ±0.4 


1.12 ±0.3 


11.8 ±3.8 


11.6 ±3.9 


32 ± 17 


33 ± 18 


10.8 



1 Missing flux determined from the 13 CO data 



1 x 10" 5 to2 x 10" 4 byl x 10~ 6 ; [CO]/[ 13 CO]: 30 to 100 
by 5; [CO]/[C 18 0]: 100 to 300 by 20). For the compari- 
son between the observed peak intensity ratios (Table 0, 
columns 3,5,7) and the predicted LVG ratios we used a 
X 2 test. To account for the absolute intensities across the 
disk of M82 we also fitted the 12 CO(J = 2 -> 1) intensity 
at each position by varying the beam filling from 0.1 to 
0.9 by 0.1. The 'best' solutions are shown for positions 3 
and 9 in Figs. || & 0- Position 3 on the western CO lobe is 
an example for a solution with low kinetic temperatures 
and high H2 densities; position 9 on the brightest MIR 
peak is representative for solutions with high kinetic tem- 
peratures and low H2 densities. 

The observed line ratios and 12 CO (J = 2 — ► 1) inten- 
sities can be modeled within the errors at all positions. 
The fit agrees very well with the data at positions where 
12 CO(J = 2 -> l)/ 12 CO(J = l -> 0) is less than 1.2. At 
position 6 (eastern MIR peak) we do not find any inter- 
section for all observed line ratios in the H2 density and 
kinetic temperature plane. For a more detailed discussion 
see Sect. LI. The best agreement with the observed line 
ratios and absolute intensities is found for a beam filling 
of 0.4. Po sition s 6 and 7 at the eastern MIR peak (Telesco 
& Gezari 1992 ) and positions 18 and 19 at the CO out- 
flow require a somewhat lower beam filling of 0.2 and 0.3 
respectively. 

The LVG parameters of the 'best-fit' across the major 



axis of M82 are shown in Figs.|]a-d. The CO abundance 
relative to H2 per velocity gradient QCO]/grad(V)) varies 
between 1 x 10~ 5 pc/kms -1 and 7x 10 _5 pc/kms _1 . As- 
suming grad(V) ~ lkms _1 pc _1 , as suggested by com- 
paring the linewidth with the linear extent of the re- 
gion, this corresponds to CO abundances in the range 
of [CO] w 10 -5 - 7 x 10~ 5 . Similar values have been 



determined in the Orion region (Blake etal. 1987) and 
were suggested by chemical models (Farquhar et al. 1994 ). 
[CO]/grad(V) increases towards the MIR peaks which in- 
dicates higher CO abundances at the active star-forming 
regions than in the more quiescent outer regions. The frac- 
tional 13 CO abundance [ 12 CO]/[ 13 CO] across M82 does 
not show any significant spatial variation. The mean value 
of all positions is 70 ± 20. A low fractional 13 CO abun- 
dance is consistent with recent radiative transfer calcula- 



tions by Mao etal. ( 2000 ) and an independent chain of 
argu ments based on CN and 13 CN measurements (Henkel 
et al. |1998|) . In contrast, the fractional C 18 abundance 
[ 12 CO]/[C 18 0] shows a trend towards higher C 18 abun- 
dances at the MIR peaks and in the outflow. While the av- 
erage [ 12 CO]/[C 18 0] ratio in the quiescent regions is about 
270, it is only about 160 at position 6, 11, 17 and 19 (see 
Fig. ||d). Note that these values suggest C 18 ab unda nces 
2-3 times higher than those used by Wild et al. ( 1992 ) for 
their LVG calculations of CO line ratios in M82. 
The kinetic temperature is well correlated with the MIR 
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emission and other tracers of high-level star formation. 
Within the prominent CO lobes, with less signs of ongo- 
ing star formation, the kinetic temperature is about 50 K. 
Towards the active star-forming regions we find two ki- 
netic temperature peaks above 150 K. These 'hot-spots' 
coincide with the location of MIR peaks (for a compari- 
son between the MIR emission and the CO distribution 



umn density at rotation level J: Nj = 93.5 



9(j+l) A J+l,J 



(1- 



see Telesco & Gezari 1992). Near the 2.2/im nucleus the 
LVG models suggest temperatures of about 75 K. Along 
the CO outflow the temperature drops with increasing dis- 
tance from the active regions. At position 17 and 19 we 
find temperatures above 100 K. At position 18 (100 pc dis- 
tance from the plane) the kinetic temperature has dropped 
to 60 K. The spatial variation of the kinetic tempera- 
ture along the major axis of M 82 is shown in Fig. || a. 
The corresponding diagram of the H2 density distribu- 
tion is shown in Fig. p| b. Solutions are found between 
n(H 2 ) — 10 2 - 7 and 10 cm~ 3 . In general, the H 2 densi- 
ties are high in regions with low kinetic temperatures and 
vice versa. The solutions for the outer CO lobes suggest 
an H2 density about n{H2) = 10 4 cm~ 3 with a tendency 
towards somewhat lower values at the very edge of the 
CO distribution (n(i?2) = 10 35 cm -3 ). These values are 
in agreement with H 2 densities calculated by Wild et al. 
( |1992| ) and Mao etal. fl2000| ). At the 'hot-spot', low H 2 



densities of n(H2) = 10 28 ~ 31 cm~ 3 are required to match 
the observed line ratios. H2 densities in the CO outflow 
are about n(H 2 ) = 10 3 cm~ 3 . 

Both the 12 CO(J = 1 -> 0) and the 12 CO(J = 2 -> 1) 
transitions are optically thick. In the cold dense regions 
we find an optical depth of Ti2 C o(j=i^o) =2 — 5 and 
Ti2 CO (j=2^i) =7 — 15. At the 'hot-spots' the derived op- 
tical depths are somewhat lower and reach unity in the 
12 CO(J = 1 — > 0) transition at the eastern MIR peak 
(position 6 & 7). For the ground transitions of the rare 
isotopes 13 CO and C 18 we find optically thin emission 
at all positions. Typical optical depths are Ti3co(J=i^o) = 
5 x 10~ 2 and t c is O (j =1 ^ ) = 5 x 10~ 3 . 



3.5. Column densities and mass distribution 

For the determination of CO and H2 column densities at 
each position we used three methods: 

- LVG: The column densities were derived from the CO 
and H2 densities, the velocity gradient and the observed 
line widths using N(CO) = 3.08 x 10 18 n(CO) gra dV v) and 

N(H 2 ) = 3.08 x 10 18 n(H 2 ) gra ^ v) , where dV is the ob- 
served line width. Therefore gri fd(v) * s an equivalent path 
length through the clouds. 

- LVGp_f (PF = partition function): The 13 CO and C 18 
column densities were calculated from the general relation 
between optical depth, excitation temperature and col- 



exp(— 4.8 \Q~ 2 v /Tex)) -1 J rdv where gj is the statistical 
weight of level J and Aj+i^j is the Einstein coefficient 
for the transition J + 1 to J. J rdv was approximated by 
/ rdv re 1.06 rdV. T ex , r, [ 13 CO] and [C 18 0] are given 
by the LVG code for each level. 13 CO and C 18 column 
densities were determined using the sum of the 6 lowest 
levels for each isotope. H2 and CO column densities were 
derived from the relative abundances of the rare isotopes 
relative to H2 and CO. 

- LTE: 13 CO and C 18 column densities were derived us- 



ing a standard LTE approach (e.g. Dickman 1978 ). As for 
the LVGpf method, CO and H2 column densities were 
derived from the abundances of the rare isotopes relative 
to CO and H^at each position. 

Column densities calculated from 13 CO and C 18 via the 
LTE method match each other with less than 5% differ- 
ence at each position. The same holds for the LVGp p col- 
umn densities calculated from T ex and r of the 13 CO and 
C 18 transition. For simplicity we therefore give in the 
following the average between the column densities cal- 
culated from 13 CO and C 18 via the LTE and LVGpf 
method. 

The spatial variations of the beam-averaged H2 column 
density across the major axis of M 82 as calculated with 
the three methods is shown in Fig. ||. The spatial distri- 
bution of the H2 column densities is in good agreement 
for all three methods. This suggests that the low J levels 
are almost thermalized. The largest difference between the 
methods is apparent at the central CO peak. While the 
LTE solutions suggest a local H2 column density maxi- 
mum of about N(H2)4" = 1 x 10 23 cm -2 , the peak is less 
prominent (N(H2)4» = 5 x 10 22 cm~ 2 ) and displaced by 
4" in the LVG and LVGpp solution (see Fig. |). 

Nevertheless, all methods clearly show that most of the 
molecular gas traced by CO is located in the outer CO 
lobes. The central 300 pc between the molecular lobes 
contain only about 20-30% of the molecular gas mass. 
Furthermore, the H2 column density distribution is 
clearly asymmetric with respect to the 2.2/im nucleus. 
We find that the centroid of mass is located about 100 pc 
south-east of the nucleus. The location of the centroid 
of mass for each method is indicated by the vertical line 
in Fig. H The highest H2 column density is found at 
the western CO lobe (position 12). Its beam-averaged 
LVG column densities are N(CO) 4 » = 2 x 10 19 cm -2 
and N(H 2 )4» = 2.3 x 10 23 cm~ 2 . The corresponding 
cloud-averaged LVG column densities are N(CO) c i ou d = 
4 x 10 19 cm- 2 and N(H 2 ) c ioud = 6 x 10 23 cnr 2 , re- 
spectively. The corresponding values for the eastern 
CO lobe (position 3) are N(CO) 4 << = 1 x 10 19 cm" 2 , 
N(H 2 ) 4 » = 1.5 x 10 23 cm" 2 , N(CO) c i oud = 2.5 x 10 19 cm" 2 
and N(H 2 )cioud = 3.8 x 10 23 cm~ 2 . For an as- 
sumed line- of-sig ht of 350 pc (for comparison with 
Mao etal. 2000 ) the mean molecular density in 



A. Weifi etal.: The state of the molecular gas in M82 



9 



Position 3 



T 12co(2-1) [K] 

30 

25 
20 
15 
10 
5 



20 



-20 



R 12 =0.98 +/- 0.28 
R 12J3 =12.1 +/" 3.7 
=44 +/- 23.7 



- R 



abundance: 12 C0/H 2 = 


7.0E-05; 12 C0/ 13 C0 = 75 


12 co/c 18 o 


= 280 


velo grad — 5.5 km s 


-l 

pc 


radiation field: 2.73 K 


beam filling 0.40 


T 12C0(1-0) = 4.1 


T 12C0(2-1) = 13. 


t 13co(i-o) = 3.2E-02 


t cibo(i-o) = 6.7E-03 


T kin = 50 K Log(nH 2 ) = 4 


cm -3 N co =1 .0E+1 9 cm" 2 




X = 1, 0.1, 0.01, 0.001 



200 

, rxi 



150 



100 



50 




3 

cc[ 



n(H 2 )] 



Fig. 6. 'Best' LVG solution at position 3 (the eastern CO lobe). The top left diagram shows the observed 12 CO(J = 
2 — ► 1) line intensities across the major axis of M 82 and the 12 CO( J = 2 — > 1) intensity as given by the LVG-model for 
position 3. The fixed LVG-parameters (abundances, velocity gradient, radiation field and beam filling) are given in the 
upper right box. The parameters below the LVG input parameters summarize the optical depth for each transition, the 
CO column density, the kinetic temperature and H2 density for the 'best fit'. The black and the grey solid line in the 
lower left diagram are the 12 CO(J = 2 -> l)/ 12 CO(J=l -f 0) (R i2 ) and 13 CO(J=l -> 0)/ 12 CO(J=1 -> 0) (Ri 2 ,i 3 ) 
line intensity ratio, respectively. The dashed lines indicate the observational errors for the corresponding line ratio. 
The C 18 0(J=1 — > 0)/ 12 CO(J=1 — > 0) (R12 is) Ihie intensity ratios are given by the dashed-dotted line. Their errors 
are indicated by the dotted line. The contours in the lower right diagram are results of the x 2 test comparing the 
predicted line ratios and the 12 CO(J = 2 — > 1) line intensity with the observed values. The star indicates the best 
solution in the kinetic temperature and H 2 density plane. 



the CO lobes is < n(H 2 ) >4" = 140 - 210cm" 3 . 
This corresponds to a volume filling factor of 
/„,4" =< n(H a ) > 4 " /n(H 2 ) « 0.01. With f aA „ = 0.4 
and a linear resolution of 65 pc we obtain characteristic 
cloud sizes of r c \ ou d = \ 65 pc / 0j 4« / f a ,4" ~ lp c - Volume 
filling factors and characteristic cloud sizes do not change 
significantly in the central star forming regions. These 
values are in good agreement with PDR models published 
by Wolfire etal. (|l990|) . 

H 2 column densities in the molecular outflow are in the 
range N(H 2 ) 4 // = 1.5 - 3.0 x 10 22 cm~ 2 . The total mass 
of the outflow is 7.2 x 1O 5 M (D = 3.9 Mpc, Sakai & 



Madore 1999) 



3.6. Conversion from I(CO) to N (H2) and total mass 

To derive the conversion factor from I(CO) to iV(H 2 ), 
we have compared LVG, LTE, and LVGpf H 2 column 
densities with the integrated 12 CO(J = 1 — > 0) intensi- 
ties at 4.2" resolution at the analyzed positions across 
the central part of M82. The variation of the conver- 
sion factor Xqo = ^V(H 2 )/I(CO) with position is shown 
in Fig. |l(^. Note that Xqq is lower than the Galac- 
tic value of 1.6 x 10 20 cm~ 2 (Kkms -1 ) -1 (Hunter etal. 



1997) at all positions and for all methods. We find that 
Xco varies across the disk of M 82 by about a fac- 
tor of 5 if one considers the LTE solutions (V"co = 
2.1 - 10.8 x 10 19 cm" 2 (Kkms -1 )" 1 ) and by a factor of 
8-9 for the LVG and \NG PF solutions {X C o = 1.3 - 

co = 1.5 - 12.2 x 



11.5 x 10 19 cm- 2 (Kkms -1 )" 1 and X, 
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Fig. 7. 'Best' LVG solution at position 9 (central MIR peak). The items and parameters are the same as in Fig[| 



10 19 cm" 2 (Kkms" 1 )" 1 ) . All methods show that the low- 
est conversion factors are associated with the central star- 
forming regions where the gas is heated by UV photons 
from the newly formed stars and cosmic-rays from SNRs. 
The CO-emitting volumes at these positions have high 
kinetic temperatures. Towards the outer molecular lobes 
with higher H2 densities and lower kinetic temperatures, 
the conversion factor rises. This is in agreement with sim- 
ple theoretical arguments that suggest that the conver- 
sion factor Xqo should be proportional to T^ n(H 2 ) 1 ^ 2 
for virialized clouds (Maloney & Black 1988). The vari- 



ation of Xqo with T,.* n(H2) 1 ^ 2 is shown in Fig. O 
The linear correlation between Xqo an d TTf- n (H 2 ) 1/2 for 
T^XHa) 172 > 0.5 is clearly visible. For T^ n(H 2 ) 1/2 < 
0.5 the scatter in the plot increases. This is in particular 
true for Xqo calculated under the assumption of LTE. 
This suggests that the gas is not close to LTE at the 'hot 
spots'. The increased scatter of A"co calculated with the 
LVG and LVGpf method might suggest that either the 
clouds are not virialized or that more appropriate models 
(like PDR models) are required to calculate the physical 
gas properties in the centre of the starburst. For a more de- 



tailed discussion see Sect. 4.1. Nevertheless, this result not 



is a function of the intrinsic gas properties which strongly 
depend on environmental effects. This implies that spatial 
variations of 12 CO(J= 1 — *■ 0) intensities can be due to 
variations of the excitation conditions of the gas rather 
than variations of column density. Similar results have 



only shows that the standard Galactic Xqo factor is not 
appropriate for a starburst system like M 82, but that Xco 



been obtained by Wild etal. (1992) using low-resolution 
CO data (see also Sect. |Q| ). Based on the analysis of Xco 
we have calculated the 'true' H2 distribution in M82 by 
interpolating the changes of Xqo from the analyzed posi- 
tions across the central CO distribution. Multiplication of 
this Xco _m ap with the integrated 12 CO(J=l — > 0) in- 
tensity distribution thus results in an H2 column density 
map. We show these maps in Fig. [l2]for Xqo derived from 
the LVG pf (top) and LTE solutions (middle) in compari- 
son with the H2 distribution one would derive assuming a 
constant, standard Galactic conversion (bottom) to illus- 
trate the importance of detailed studies of Xqo to derive 
H2 column density distributions. The H2 column density 
maps in Fig. |l2| (top & middle) indicate that the central 
star-forming region is surrounded by a double-lobed dis- 
tribution of molecular gas, while H2 seems to be depleted 
in the central starburst region itself (see also Fig. ||). 
The total H2 mass of the region shown in Fig. is 
2.3 x 10 8 M Q for the ING PF and LVG and 2.7 x loHl 
for the LTE solution at a distance of D = 3.9 Mpc (Sakai & 
Madore 1999). The corresponding values at D = 3.25 Mpc 
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Fig. 8. LVG solutions for positions 1 to 16. Top: locations of the analyzed positions. The radii of the circles indicate the 
spatial resolution for which the line ratios have been determined, a) to d): spatial variations of the kinetic temperature, 
the H2 density, the CO abundance per velocity gradient and the fractional C 18 abundance across the major axis 
of M82. The error bars in a) and b) correspond to the parameter range of kinetic temperatures and H2 densities for 
which the LVG line ratios and the 12 CO(J= 2 — > 1) intensity is consistent with the observations within the errors. 
This corresponds to the area within the \ 2 — 1 contour shown for position 3 and 9 in Figs^ and ^ . The error bars in 
c) and d) correspond to the range for [CO]/grad(V) and fractional C 18 abundance where the x 2 of the corresponding 
fit is 50 times higher than for the 'best' solution. 



(Tammann & Sandage 1968) are 1.6 and 1.9 x 1O 8 M0, J^.l. Comparison with other radiative transfer analyses 



respectively. These values are in good agreement with 
estimates from 450 /im dust continuum measurements 
(Smi th et a l. |l99l[ ) and from C 18 0(2^1) intensities (Wild 



etal. 1992). Therefore, the total molecular mass is 3 times 



lower than the mass one would derive using the standard 
Galactic conversion factor of 1.6 x 10 cm" - 2 (Kkms" 1 )- 1 
(4.9 x 10 8 M D = 3.25 Mpc; 7.1 x 10 8 M Q D = 3.9 Mpc). 



4. Discussion 



Analyses of the physical conditions of the molecular gas in 
M82 have been published by Tilanus etal. fll99l|), Wild 



etal. (1992), Giistcn etal. ( |1993[ ) and more recently by 
Mao etal. ( |2000D and Petitpas & Wilson ( ^000[ ) using 
single-dish CO data and other tracers of the molecular 
gas. The kinetic gas temperature of the CO-emitting gas 
phase derived in these studies are of order Tfe in = 30 — 
70 K. H2 densities range between n(H2) = 10 3 — 10 4 cm~ 3 . 
Thus our solutions at the CO lobes and the outer parts 
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of the CO distribution (n(H2) = 10 35 ~ - cm -, ±& 
45 — 60 K) are consistent with previous studies. 
The situation is different for our LVG solutions within the 
starburst region. Kinetic gas temperatures above 150 K 
are clearly inconsistent with results published in litera- 
ture so far. In the most recent analysis of the excitation 
conditions of the molecular gas using mm and sub-mm 
CO emission lines (up to J 



(2000) 



7 - 6), Mao etal 
suggested kinetic gas temperatures as high as 130 K. How- 
ever, they rejected their LVG solution because of intrin- 
sic inconsistencies regarding almost equal area and vol- 
ume filling factors derived from the one-component LVG 
model. This leads to characteristic cloud sizes of about 
150 pc which is inconsistent with high-resolution studies 
of the CO d istribution in M 82 (Shen & Lo |l995| , Neininger 
etal. 1998 and this work). Characteristic cloud sizes de- 
rived from our analysis are only w 1 pc, hence more real- 
istic. This difference results mainly from a very low area 
filling facto r of only f a<2 2" = 0.04 - 0.07 found by Mao 
etal. (2000). From the CO morphology (see Fig. P and 
fa, 4" = 0.4 — 0.2 we would expect area filling factors of 
fa, 22" = 0.2 — 0.15 at 22" resolution. The reason for this 
discrepancy remains unclear. But obviously the assump- 
tion of an isothermal gas phase used in the LVG model 
is more reasonable for our high spatial resolution study 
than for the low-resolution data used by Mao et al. ( 200(3 ) . 
Even though our LVG analysis does not lead to internal 
inconsistencies we also find that the observed line ratios 
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Fig. 9. Beam averaged H2 column densities across the ma- 
jor axis of M82. Offsets are given relative to the centre of 
M82 (2.2/im peak). The thick solid line corresponds to 
the H2 column densities derived from the LVG solutions. 
The dashed and the dotted lines give the corresponding 
values for the LTE approximation and the solution calcu- 
lated from the partition functions (LVGp F ), respectively. 
The vertical lines indicate the location of the centroid of 
mass of the three H2 distributions. Note that the loca- 
tion of the centroid of mass is identical for the LVG and 
LVGpf solution 

. The labeling of the right axes gives the H2 mass per beam in 
units of solar masses. 
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Fig. 10. Variation of Vco across the analyzed positions 
in M82. Offsets are given relative to the centre of M82 
(2.2/im peak). The H2 column densities were derived us- 
ing the LTE, the LVG and LVG PF solutions. The dotted 
line at Xco — 1-6 x 10 20 cm -2 (Kkms -1 ) -1 corresponds 
to the standard Galactic conversion factor. 
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Fig. 11. V co as derived from the LTE, LVG and LVG PF 
solutions versus T _1 n^Bk) 1 ^ 2 . The plot shows a linear 
correlation between the conversion factor and the as- 
sociated H2 densities and kinetic temperatures beyond 
T- 1 !!^) 172 « 0.5. 



are difficult to reproduce with the one-component LVG 
model at the 'hot spots'. This is in particular true for 
positions 6, 7 & 11 (eastern MIR peak and expanding su- 
perbubble) where no intersection of all observed line ratios 
(disregarding the errors of the observations) exists within 
the calculated parameter space. At these positions more 
sophisticated radiative transfer models like PDR models 
probably lead to more consistent results. However, a com- 
parison between LVG and PDR models in M 82 shows that 
constraints on H2 densities and beam-averaged column 
densities are very similar for both methods (Mao et al. 
|2000D . Giisten etal. ( p.993fl and Mao etal. Q2000D con- 
cluded that in order to explain the observed line ratios, 
a two-component model of the molecular gas in M82 is 
needed. However, densities of n(H2) = 10 3 7 cm -3 for the 



high density and n(H.2 ) = 10 cm 3 for the low density 



component (Mao etal. |2000 ) are much lower than values 
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Fig. 12. H2 column density maps calculated from the in- 
tegrated 12 CO(J = 1 — > 0) intensity distribution at 4.2" 
resolution. The H2 maps were derived using the Xqq dis- 
tribution across the major axis of M 82 as calculated from 
the LVGpf solutions (top), from the LTE solutions (mid- 
dle) and from the standard Galactic Xqq, conversion fac- 
tor (1.6 x 10 20 cm- 2 (Kkms -1 )- 1 ) (bottom). The total 
H2 mass is denoted in the upper right corner of each plot. 
Contour levels correspond to 1.4, 2, 3, 4, 6, 8, 10, 12, 
14 and 16 x 10 22 cm" 2 (top), 3, 4, 6, 8, 10, 12, 14 and 
16 x 10 22 cm~ 2 (middle) and 7, 10, 14, 18, 22, 26 and 



30 x 10 22 1 



(bottom). 



given by Gusten etal. ( |1993| ) (10 3 - cm -3 , 10 5 -°cm" 3 ). In 
addition the high density component in the models by 
Gusten etal. is associated with a high kinetic tempera- 
ture while it is cooler than is low density component in 
the models by Mao et al. ( 2000| ). Our high-resolution LVG 
results confirm the PDR calculations of Mao et al. ( 20001 ). 
In particular, their conclusion that the bulk of CO emis- 
sion in the core of M 82 arises from a warm, low-density 
interclump medium is consistent with our findings. 
Further support for a highly excited CO component to- 
wards the active regions in M 82 comes from the mor- 



phology of the high- J CO lines observed by Mao et al. 
( 2000D . They find that the spatial separation of the CO 
lobes decreases with increasing J. While the spatial sepa- 
ration of the outer CO lobes in the 12 CO( J — 1 — > 0) and 
12 CO(J = 2 — * 1) transition is about 26", it decreases to 
only 15" in the 12 CO(J = 7 — > 6) transition. This distance 
is in good agreement with the spatial separation of the 
kinetic temperature peaks that we find in our LVG solu- 
tions (Fig|a). 

To further test the reliability of the modeled kinetic tem- 
perature and H2 density distribution across M 82 we cal- 
culated line ratios for the high- J transitions of 12 CO 
and 13 CO at 22" resolution and compared o ur pre diction 
with the line ratios published by Mao etal. ( 200C ). Note 
that the spatial smoothing of our high-resolution one- 
component excitation model leads to a multi-component 
model at lower resolution because it encompasses the indi- 
vidual solutions (weighted with a Gaussian of 22" width) 
at all positions. The predicted line ratios from a single CO 
isotope match the observations extremely well. For differ- 
ent CO isotopes the predicted line ratios are lower than 
suggested by the observations, but consistent within the 
observational errors. The predicted and observed line ra- 
tios are summarized in Table |3|. 

An independent chain of arguments favoring a tempera- 
ture gradient towards the central starburst region comes 
from the different separation of the eastern and western 
'hot spots' in MIR and FIR observations. Hughes et al. 



(1994) stated that the larger separation of the peaks at 



450 /im reflects the radial temperature gradient that must 
exist within the torus if the dust is heated by the central 
starburst population. 

4.2. The state of the molecular gas in M 82 's centre 

Our study does not imply that the central part of M 82 
does not contain dense molecular cloud cores as traced by 
other molecular species lik e HC N (see e.g. Seaquist et al. 
Brouillet & Schilke |l993|). Most of the CO emis- 



1985 



sion from the core of M82, however, arises from a warm, 
low density interclump medium. This conclusion was al- 
ready reached by Mao etal. ( 2000| ) using PDR models. 
The existence of such a gas component, which we can be 
confident of, raises the question if it can survive under the 
extreme conditions in the starburst centre or whether it is 
indicative of cloud evaporation and thus of a destruction 
of molecular clouds. From the observed UV field strength 
one would expect that the molecular clouds with densities 
less than 10 4 cm~ 3 cannot exist in M82's centre (Brouil- 



let & Schilke 1993). In addition, such a gas component is 
barely dense enough to avoid tidal disruption (Mao et al. 
|2000 ) which would support a rather diffuse distribution of 
the gas. Such a scenario would explain the high [CI] /[CO] 
abundance ratios observed in M 82 (Schilke et al 



1993 



White et al. 1994 ) and the depletion of H2 in the central 
300 pc of M82. As a consequence, this would imply that 
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Table 3. Line ratios for the high-J transition of 12 CO and 
13 C0 at the centre and the eastern and western CO lobe 
of M82. The line ratios are given for a spatial resolution 
of 22". The first row of each line ratio gives the value 
calculated from the LVG model. The second and third row 
are the observed peak and integrated line ratios adopted 
from Mao etal. ( |2000[ ). 



CO line ratio 




east 


centre 


west 




12 CO(7->6)/ 12 CO(4- 


-3) 


0.36 


0.25 


0.31 


model 






0.32 


0.29 


0.37 


peak 






0.28 


0.36 


0.33 


integral 


12 CO(4^3)/ 13 CO(3- 


.2) 


7.07 


8.02 


5.81 


model 






8.50 


10.1 


9.07 


peak 






9.09 


6.88 


7.40 


integral 


13 CO(2->l)/ 13 CO(l- 




1.92 


1.67 


1.87 


model 






1.53 


1.63 


1.93 


peak 






1.47 


1.44 


1.17 


integral 


12 CO(2-^l)/ 13 CO(2- 


-1) 


7.82 


9.27 


7.48 


model 






9.04 


13.7 


9.19 


peak 






11.4 


14.2 


11.7 


integral 



the clouds are not virialized. This could explain why the 
linear correlation between Tj^ n(H2) 1 ^ 2 and Xco is not 
valid for n(H 2 ) 1/2 < 0.5, hence for the warm and low- 
density gas in the star burst centre (see Fig. |ll|). Therefore, 
there seems to be observational evidence that molecular 
clouds are partly disrupted and dissociated in the star- 
burst centre. More detailed numerical analysis of evap- 
oration time scales and related questions are needed to 
settle this question. 



4-3. Xqo 

The suitability of a global Galactic factor Xco to con- 
vert 12 CO(J = 1 — > 0) intensities to H2 column densi- 
ties has been discussed by many authors (e.g. Young & 

1984 Hunter etal. [L997D. As 



Scoville 1982, Bloemen etal 
seen, theoretical studies of Xqo 



showed that it is sensi- 
tive to the kinetic temperature of the emitting gas and 
that the conversion factor should be lower for starburst 
galaxi es like M82 than for the Milky Way (Maloney & 
Black 1988). Investigations of Xqo m M82 confirmed 



this prediction: Wild etal. (1992) used the optically thin 
C 18 0(2— >1) transition to derive H2 column densities and 
hence Xqo along the major axis of M82. They found 
Xco < 1 x 10 20 cm -2 (Kkms -1 ) -1 and variations by a 
factor of 2 along the major axis. Similar results were ob- 
tained by Smith et al. fll991| ) using the 450 /xm continuum 
emission from dust grains to derive H2 column densities. 
Even though both studies suggest a low conversion factor, 
its variation across the major axis shows significant differ- 
ences. While Smith et al. ( |l99l| ) found that Xco decreases 



from east to west, with no changes in the central star- 
forming regions, the results by Wild etal. ( 1992j ) suggest 
very low conversion factors near the eastern MIR peak 
and an increasing Xco towards the western molecular 
lobe. Therefore our variations of Xco i n general support 
the results by Wild etal. (1992). Nevertheless our conver- 



sion factors are slightly lower than those inferred by Wild 
etal. ( 1992| ), and the location of the western Xco maxi- 
mum is displaced by w 7". The discrepancies between the 
Xco factors derived from the molecular lines (Wild et al. 



1992 and this work) and the estimates from the dust emis- 
sion might resul t from the single-temperature model used 
by Smith etal. ( 1991 ). In particular in the central region, 
which shows strong MIR emission from heated dust (Tele- 
sco & Gezari 1992 ), this might lead to an overestimate of 
the H2 column density and thus to an overestimate of 
Xco • Furthermore, the different morphology of the 450 
fim map published by Hughes et al. ( |1994| ) raises doubts 
on the reliability of the 450 /im intensities used by Smith 
etal. (1991) for their calculation. From this we conclude 
that the Xco factor in M82 is not only lower than the 
standard Galactic conversion factor, but that in addition 
Xco in the central 300 pc is at least 3 times lower than 
in the molecular lobes. A similar gradient for the conver- 
sion factor has been found in the Milky Way towards the 



Galactic Centre (e.g. Blitz etal. 1985, Sodroski etal. 1994, 
Dahmen etal. 1998). Furthermore, our analysis suggests 
that the variations of Xco ar e mainly caused by variations 
of the kinetic temperature of the CO-emitting volume due 
to environmental effects while abundance variations play 
a minor role. 



5. Conclusions 

We have observed the 12 CO(J = 2 -> 1) and 
C 18 0(J = 1 — ► 0) emission lines in the starburst 
galaxy M82 with high spatial resolution using the 
Plateau de Bure interferometer. Our main conclusions 



1) The overall morphology and kinematics for both 
transitions are similar to those of 12 CO(J = 1 — > 0) 
and 13 CO(J = 1 -> 0) pu blished by Shen & Lo fll995| ) 
and Neininger etal. (1998). The dynamical centre of the 
molecular gas coincides with the 2.2 /jm nucleus while 
the centroid of the molecular mass is located 100 pc 
west of M 82's centre. South of the expanding molecular 
superbubble (Weifi etal. [1999 ) an outflow of molecular 
gas with a total mass of 7.2 x 10 5 M is detected. 

2) The 12 CO(J = 2 -> l)/ 12 CO(J=l -> 0) line intensity 
ratios are lower (< 1.4) than previously reported. Thus, 
CO line ratios in M82 are not outstanding, but compa- 
rable . to values found in other starburst galaxies like 
NGC253. Line ratios vary across the disk of M82. Near 
the MIR peaks, the 12 CO(J = 2 -> l)/ 12 CO(J=l -► 0) 
ratios are high; in the outer parts, that are less affected 
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by the starburst, this ratio drops to unity. 
3) An LVG excitation analysis of the CO lines suggests 
that the excitation conditions of the molecular gas are 
strongly influenced by environmental effects. In the 
outer parts of the CO distribution we find H2 densities 
of n(H2) = 10 35_42 cm~ 3 and kinetic temperatures of 
Tfem = 45 — 60K. Towards the star-forming regions, 
indicated by strong MIR emission, the kinetic tempera- 
tures raise above 150 K. The hot gas is associated with 



low H2 densities of only n(H2) 



10 



2.8-3.1, 



Area 



filling factors of f a = 0.2 — 0.4 and volume filling factors 
of /„ = 0.001 — 0.02 indicate that the gas is organized 
in small clumps with a typical size of r c i ouc i w 1 pc. 
[ 12 CO]/[ 13 CO] abundance ratios are about 70 with- 
out significant spatial variations across the galaxy. In 
contrast, [ 12 CO]/[C 18 0] abundance ratios in the outer 
parts of M82 are comparable to those found at the 
Galactic Centre ([ 12 CO]/[C 18 0] = 270) but decrease 
to only [ 12 CO]/[C 18 0] = 160 at the star-forming re- 
gions. Beam-averaged H2 column densities range from 
N(H 2 ) 4 " = 2.4 x 10 22 cm~ 2 near the MIR peaks to 
N(H 2 ) 4 » = 2.3 x 10 23 cm~ 2 at the western CO lobe. The 
H2 distribution has a double-peak morphology which 
surrounds the central starburst region. The central 300 pc 
are depleted in H2. Thus the H2 distribution differs from 
the CO distribution. This result even holds when the H2 
column densities are calculated under the assumption of 
LTE conditions. The total molecular mass is 2.3 x 10 s M Q . 
4) The conversion factor from /(CO) to N(H 2 ) (X C o) 
depends on the excitation conditions of the CO-emitting 
volume. Even in regions which are less affected by the 
starburst A"co is about 3 times less than the standard 
Galactic value. From the LVG analysis we find that 
Xco ~ n(H2) 1 ^ 2 . Therefore Xqo is lower in the 

central star-forming regions than in the outer molecular 
lobes. 

Appendix A: Short Spacing Correction 

A major difficulty related to flux determinations from in- 
terferometer maps is the missing flux problem. In brief it 
arises from the lack of coverage at low spatial frequencies 
in interferometric observations. This leads to an insen- 
sitivity to emission more extended than A/S^f, where 
^min i s tne length of the shortest projected baseline. The 
recovered flux therefore only represents clumpy parts of 
a brightness distribution which leads to an underestimate 
of the total flux of a source. For a more detailed descrip- 



tion see e.g. Wilner & Welch 1994 and Heifer & Blitz 



1995. The only way to overcome this problem is to derive 



the visibility at the origin of the ui>-plane (which repre- 
sents the integrated flux of a brightness distribution) and 
the low spatial frequencies from single-dish measurements 
and combine them with the interferometric observations. 
Methods for the combination of single-dish and interfer- 
ometer data have been described e.g. by Vogel et al. ( 1984 ) 



and Herbstmeier etal. ( 1996| ). Both methods generate the 
central visibilities from the single-dish data, as if they 
were measured with the interferometer. The combined vis- 
ibility set is then processed in the standard interferometer 
reduction procedure. In this method the relative weight 
between single-dish and interferometer visibilities has a 
strong effect on the resulting brightness distribution. Fur- 
thermore the dirty image calculated from the combined 
visibility set still needs to be deconvolved using CLEAN or 
other deconvolution algorithms. In particular the CLEAN 
deconvolution algorithm is problematic for short-spacing 
corrected images because it fails on extended structures. 
We therefore used a different approach to combine the 
PdBI and 30 m telescope data. 

The basic idea behind our method is that the missing flux 
problem only arises from an incorrect interpolation of the 
visibilities in the central part of the wy-plane (V u 2 + v 2 < 
A/S^f). Therefore the missing flux problem in a finally 
reduced interferometer map can be solved by replacing the 
questionable part of the Mi^-plane by the values calculated 
from a single-dish map with identical extent, grid and flux 
units. This procedure avoids additional CLEANing on the 
combined data set and the choice of different weightings 
between interferometer and single-dish data. The require- 
ment for the single-dish data is the same as described 
by Vogel etal. ( |1984| ). For the combination we regridded 
the single-dish data cube to the same spatial and veloc- 
ity grid as the interferometer data. We then converted the 
flux units from T m \, to Jy/pixel using the Ray leigh- Jeans 
approximation, S/beam = jsT m b and beam — 1.133 / 2 
where f s is the FWHM of the 30 m telescope beam in units 
of the pixel size. The flux units of the cleaned and primary 
beam corrected interferometer cube were also converted to 
Jy/pixel with beam = 1. 133 fi ; maj-fi,min, where fi, ma j and 
fi,min are the FWHM of the major and minor axis of the 
clean beam in units of the pixel size. Furthermore we gen- 
erated a model for the 30 m telescope main beam and the 
interferometer clean beam. The 30 m beam was assumed 
to be represented by a circular Gaussian with FWHM = 
f s . The interferometric beam was described as a Gaus- 
sian with major and minor axis fi, m aj and fi >m i n . The 
normalization of both Gaussians was such that the am- 
plitude of the visibility at the origin of the ui>-plane was 
1. We then transformed both data cubes and the model 
beams to the wt>-domain using an FFT algorithm. The real 
and imaginary parts of the single-dish data were divided 
by the amplitudes of the model 30 m beam to deconvolve 
the single-dish visibilities from the 30 m telescope beam. 
In order to match the interferometer data the result was 
then multiplied by the amplitudes of the clean beam. At 
this stage of the combination the real and imaginary parts 
in the single-dish and interferometer data are comparable 
and the central interferometer pixels can be replaced by 
the single-dish values. The part of the wu-domain to be 
replaced by the single-dish values in general depends on 
the spacing of the shortest baseline and on the effective di- 
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Fig. 13. Flow chart of the short-spacing correction. 

ameter of the single-dish telescope. For data sets with no 
overlap in the wif-domain we replaced the part that corre- 
sponds to the effective diameter of the 30 m telescope (am- 
plitude of the visibilities for the 30 m beam model > 0.5). 
Otherwise we selected the part smaller than the short- 
est interferometer baseline. Finally the combined real and 
imaginary parts were transformed back to the image do- 
main and the flux units were converted to Kelvin consid- 
ering that the combined beam is equal to the clean beam. 
Note that there are no free parameters in the combination 
of the data sets except for the choice which part of the uv- 
plane is replaced by the single-dish data. The methods re- 
quire the knowledge of the single-dish beam pattern and 
the cle an beam only. A flow chart for the Short-Spacing 
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